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The quantum Fourier transform (QFT) is a powerful tool in quantum computing. The 
main ingredients of QFT are formed by the Walsh-Hadamard transform H and phase shifts 
P(-), both of which are 2x2 unitary matrices as operators on the two-dimensional 1-qubit 



CN , space. In this paper, we show that H and P(-) suffice to generate the unitary group U{2) and, 



consequently, through controUed-f/ operations and their concatenations, the entire unitary 



I group U{2'^) on n-qubits can be generated. Since any quantum computing algorithm in an 
^ ■ n-qubit quantum computer is based on operations by matrices in f/(2"), in this sense we 
■ have the universality of the QFT. 
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1 Introduction 

The quantum Fourier transform (QFT) on the additive group of integers modulo 2™ is 
defined by 

2™--! 

^2™(|a))= ^e(2-^'^^)/2'"|l/), for a G {0, 1, 2, . . . , 2" - 1}. (1) 

y=o 

QFT plays a significant role in the development of the quantum computer (QC). One may 
note, for example, that the potentially powerful integer factoring algorithm by P. Shor relies 
critically on the QFT for the detection of periodicity springing from the prime factors. 
We can further analyze (|I|) as follows. First, write 



a = ai2'"~^ + 022™"^ H h am-i2^ + 0^2° = (aiOa . . . 



and 



y = yi2™-i + 2/22™-' + ■ ■ ■ + y^_i2^ + ?/„2° = {ym ...y„ 
Then it is well known that 



RHS of (1)= ^e(2'^^'^^/2'")|2/i...U 

^2TTi(0.am)yi |„, \ 2m{0.am-iam)y2 



_ "y ^ g27ri(0.a™)j/i |^^^g27ri(0.am_ia™)2y2 1^^^ . . . g27ri(0.aia2...am)ym 1^ ^ 

y=o 

In the above factorization (or "untangling"), each factor is of the form 

|0) + e-|l). (3) 
Such a state can be produced in two steps P, pp. 340-341]: First, apply the transformation 



V2 



1 1 
1 -1 



(4) 



where H is known as the Walsh-Hadamard transform, to the state |0): 



^^|0) = ^(|0) + |1)). 



(5) 



Next, apply the phase shift operator 



Pioo) 



1 

e'"^ 



(6) 
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to yielding 



Pico)[Hm = —i\0) + e^-\l)). (7) 

The RHS of (|^) is (apart from a normalization coefficient). Therefore, we see that the 
constituents of the QFT are H and P{uj). From the quantum optics point of view, H is 
realized by a half-silvered mirror (beam splitter) and P{uj) represents a phase shifter, as in 
a standard Mach-Zehnder interferometer ([^ Q). 

First, we wish to emphasize that the QFT strictly by itself is not universal in quantum 
computing; see Remark |^ below. Thus, the question becomes whether the two constituents 
H an P(-) of QFT are universal or not. The conjecture we want to pose here is the following: 

[Q] "Any QC algorithm can be represented as a composition of 

Walsh-Hadamard transforms and associated conditional phase shifts." (8) 

The implication of (H) is that the realization of any QC algorithm translates into a com- 
bination of elementary quantum interferometric operations, i.e., single particle beam splitter 
(Walsh-Hadamard transform) followed by a conditional phase shift. Any QC algorithm can 
thus be formulated, or reformulated, in terms of elementary multiparticle quantum interfero- 
metric operations. The unique universal fundamental properties of QC concerning quantum 
superposition, entanglement and interference are all explicitly represented in terms of quan- 
tum multiparticle interferometry (QMI). 

QMI practically is not to be taken as a proposed embodiment of a QC any more than the 
Turing machine is to be taken as a literal construction in classical computing. Rather, Ekert 
1^ has suggested its equivalence to QC in the sense of its universality, meaning that QMI 
could be viewed as the closest QC analogue of the classical Turing machine (through the 
universality theorem established in this paper). This concept and viewpoint should provide 
physical insights into the operational aspects and can facilitate efficient design of a universal 
QC. 

2 Mathematical Proof of the Universahty of H and P(') 

Our answer to [Q] is affirmative. We now proceed to provide the mathematical justifica- 
tions below. 

As usual, we let U{n) to denote the unitary group on n-dimensional space. By abuse of 
notation, we regard U (n) the same as the multiplicative group of all n x n unitary matrices. 
SO{n) denotes the orthogonal group on n-dimensional spaces or, equally, the multiplicative 
group of all n X n orthogonal matrices. We also define the maximal torus T{n) in U (n) as 

T(n) = {diag(e'"S . . . ,e'"") |c^i,c^2,... ,c^neR}, 

i.e., T[n) consists of all n x n diagonal matrices whose diagonal entries are complex numbers 
of unit magnitude. T(n) is a subgroup of the multiplicative group U{n). 
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Let ^ be a collection of n x n unitary matrices. In this paper, we will use Qni-A.) to 
denote the unitary subgroup ofU{n) generated by A, i.e., 

Gn{A) = [^{Gq I Gq, is a subgroup of U{n),A C Ga}- 

a 

We will write Qn{A) simply as Q{A) if the value of n is clear from the context. 
We begin with n = 2. 

Lemma 1 Lemma 4.1]). We have U{2) = g{SO{2),T{2)), i.e., U{2) is generated by 
50(2) and T{2); more precisely, for every A & U{2), we have 



A 



,i5 Q 



e'^ 
for some a, /3, (5, G M. 



gia/2 Q 

e-^"/2 



cos u sm u 
— sin oj cos oj 







Lemma 2. T{2) C g{H,P{-)). 

Proof. We first note that the NOT-gate 

X 



1 

1 



can be obtained as 

X = HP{-tt)H. 
Therefore X G Q{H,P{-)). From this, we have 



' 


1 " 




" 1 " 




' 


1 " 




' 1 " 




- Q 


1 







e*"i 




1 







e*'^2 




e'"^^ 



XP{uJi)XP{iU2) 

for any given uji,uj2 £ I^- Therefore Q{H,P{-)) contains the maximal torus T(2). 

Lemma 3. S0{2) C g{H,P{-)). 
Proof. For each rotation matrix 



R{uj) 



cos LU sm LU 

— sin uj cos uj 



we easily verify that 



R(uj) = P (II) HP{uj)XP{-uj)HP (|- 



□ 



(9) 
(10) 

(11) 
□ 



(12) 
□ 



Theorem 4. g{H,P{-)) = U{2). 
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Proof. This follows immediately from Lemmas 



□ 



At this point, it should already be clear from the results in that f/(2") can be gen- 
erated through controlled-U [2) gates, for any n = 1,2, . . . . To make this paper sufficiently 
self-contained, however, let us give the following concise, rigorous treatment as to how to 
construct any V G [/(2") from a serial connection of a collection of unitary matrices Vij, 
where each Vij is a (generalized) controUed-f/ (2) gate. The precise statement is given below. 



Theorem 5. Let V E t/(2"). Then 



V 



2"-l i-1 
i=l j=0 



(13) 



for a collection of matrices Vij G U{2"') such that 



V,f. S., 



Sij is the identity transformation, 



Sij = span{\m) | m G {0, 1, . . . , 2"- — 1} , m ^ i, m ^ j}, 
< J < i < 2" - 1. 



(14) 



In other words, each V G [/(2") is a product of (generalized) controlled-U (2) unitary matrices 
Vij, which acts nontrivially only on S^- = span{\i), \ 

Proof. We first quote the following fact ||5|, ||: For any V G U{2'^), there exists a collection 
of unitary matrices Tjj, < j < i < 2" — 1, and a D G T(2"') such that 



^1 i-l 



V 



n n^M 

i=l jr=0 



D, 



(15) 



where Tjj G 5*0 (2") C f/(2"') is a rotation involving \i) and \j) and satisfying (0). For the 



benefit of the reader and for the sake of self-containedness, we include a direct proof of (15) 
in the Appendix, condensed from 



Now we can break up D into 
D = 



di 



\ 



\ 



(16) 



where 



/ d, 




di 



(17) 



5 



and 



/ 1 



di 



(18) 



V 



1 / 



for i = 2,3, .. . ,2" — 1. It is easy to see that Di acts trivially except on |0) and and 
the other Dj's act non-trivially only on In addition, Dj's commute with each other, and 
each Di commutes with VO < / < A; < i as well. Thus, 

V = T2n_i 2"-2 • • • ^2"-l,0^2"-2,2"-3 • • • 72"-2,0 • • • 2"-l 

= 72n_i^2"-2^2"— 2,2"-3 • • • ^2"-l,0^2"-l 
^2"-2,2"— 3 • • • ^2"-2,0-O2"-2 



> 2" — 1 strings of products 



^2,1^2,0-02 

For < J < z < 2" - 1, define 



(19) 



T- ■ 
TijDi 



if J =0. 



Therefore we have reached 



V 



2"-l i-1 
i=l j=0 



where each Vij is a unitary matrix which acts nontrivially only on the states \i) and \j) 
satisfying (0). □ 



Remark 1. (1) In Barenco et al. [Q, p. 3465, right column, line 34], the equation there 
corresponds to our equation (|T5|) here. However, a summation sign ^ is used instead 

of the product sign Y\ (which is actually a double product H 11 i'^ dH)) which, of 

j i 

course, is a misprint. 

(2) The factoring of D in ([l^) into the product of Di,D2,... and -D2"-i in the form of 
( p!7D and (|18|) is peculiar in the sense that Di is chosen differently from the other -Dj's, 
i ^ 1. It must be done this way (but no further mathematical explanations were given 
in Ijl]]). The reason for this is that there are 2" — 1 strings of products as indicated in 
(|19]). Therefore D must be factorized to have 2" — 1 factors Di, D2, . . . , D2"-i, in the 
unique way of (|1^ and (p!8| ) in order to satisfy (|T^). □ 
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Remark 2. Now it can be readily seen that the QFT itself is not universal in the sense that 
f/(2") is not generated by (cf. (0), with m = n therein) or (generalized) controlled- JF2m 
(where m < n) operations. First, check n = 1: we see that J^2^ = ^2 is actually the 
Walsh-Hadamard transform H (apart from the normalization factor 1/a/2). Therefore, the 
phase shifts P{oj) in (|^) cannot be generated by J-2 because P{uj) has eignevalues 1 and e*"^ 
while H has eigenvalues 1 and —1. For a general positive integer n, the range of J-'2" or of 
controlled- ^-"2™ , m < n, consists at most of linear combinations of states of the form 

g2^i[{0.a„)yi + (0.a„_ia„)y2+...+(0.ai...a„)y„]|^^ _ _ ^^^^^ ^,^y, ^ j = 1,2, . . . ,n. 

The phases of such states are not even dense with respect to all possible phases e^'^*^, < 
e <2tt. □ 



3 Remarks on Circuits 

The decomposition (^) is a mathematical rendering of statement [Q] and answers the 
conjecture affirmatively. In this section, let us further elaborate on the circuit design aspects, 
based on the work in Barenco et al. [|l], §VIII] and Q]. 

) satisfies ([T^) and thus Vij acts nontrivially only on the states \i) 
and Denote the restriction of Vij to the 2-dimensional subspace = span{|i), |j)} by 

Vij. Then Vij G U{2). As pointed out in [|l], p. 3465], each Vij is not a standard A„_i(Vij) (in 
the notation of [||, p. 3458]) gate in the sense that the controls are states rather than bits. 

Nevertheless, using Proposition ^ below, Barenco et al. §VIII] point out how to rear- 
range basis states with a "gray code connecting state \i) to state such that Vij becomes 
unitarily equivalent to A„_i(Vij). In this sense, Vij are generalized controlled- gates. 

Proposition 6. The symmetric group S2" of permutations on the symbols 0, 1, 2, . . . , 2" — 1 
IS generated by the 2-cycle (2" - 2, 2" - 1) and the -cycle (0, 1, 2, . . . , 2" - 1). 

Proof. This is a basic fact which can be found in most basic algebra or group theory books. 

□ 

Incidentally, we note that the 2-cycle (2" — 2, 2" — 1) is a permutation between the states 
I 1 1 ■ — 1 Q ,) and I 1 1^. . 1 ) and thus can be realized by the controlled-NOT gate with the 

n bits n bits 

nth qubit as the target bit and the first (ra — 1) bits as the control bits as shown in Fig. 1. 

On the other hand, the 2^-cycle (0, 1, 2, . . . , 2" — 1) makes the rotation of the states 

|0) ^ |1) ^ ^ |2" - 2) ^ |2" - 1) ^ |0), i.e., the |a;) ^\x + l mod 2") operation. This 

can be implemented by the circuit as shown in Fig. 2. 

Therefore, any permutation of the basis states a; = 0, 1, 2, . . . , 2" — 1, can be realized 
by finitely many controlled-NOT operations consisting of circuits as shown in Figs 1 and 2. 

Thus, each factor Vij in (^) can be realized by the circuit as shown in Fig. 3. 



Each factor Vij in (13 
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X 



Fig. 1 The n-bit controlled-NOT gate A„_i(X), where X is given by (P). This 
gate implements the two cycle (2" — 2, 2" — 1) in Proposition ^. 
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Fig. 2 This circuit implements the operation |^) — » |x + 1 mod 2") or, equivalently, 
the 2"-cycle (0, 1, 2, . . . , 2" - 1) in Proposition §. Note that the bit |1) at the 
bottom of the figure is the "scratch bit" which is sometimes omitted in circuit 
drawing. All the gates in this circuit are controlled-NOT gates. 
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Fig. 3 The unitary matrix Vij in ( p!B| ) as a controlled- Vi^ gate where V^j € U{2). 
The operations (z, 2" — 1) and (j, 2" — 2) in the two boxes are cyclic permutations 
(which can be realized by concatenations of circuits in Figs. 1 and 2). 

By concatenating together all the blocks Vij as shown in Fig. 3 according to the factor- 
ization (p!3D , we have constructed all V G [7(2") with controlled- V^j gates according to ([131) . 
Each Vij G U{2) is then further formed from concatenations of the gates H, P{uj) G U{2) by 
Theorem ^. It is in this sense that we have the universality of the Walsh-Hadamard gate 
H and the phase shift gate P(-) and, consequently, that of the quantum Fourier transform 
with the affirmative answer to question [Q] in (|Tot). 



Appendix: Decomposition Procedure of General Finite 
Dimensional Unitary Transformations into a Product of 
Plane Unitary Transformations 

First, we define a special type of unitary transformations Tpq(0,cr) G U{n) by 
TrJ(J), a) = [tij]nxn, l<P,q<n, q, 



where 



1, 
COS( 

.0' 

e'" sin I 



i = j,i 7^ P,i 7^ q, 

i = j= p or i=j = q, 

i ^ 3,^ ^ PJ ^ q ^T^di ^ qj ^ p, 

i = p and j = q, 

i = q and j = p; 



9 



I.e., 



P 



Tpq{(l>, cr) = 



10 
1 

1 



COS0 — e ^^'sin^ 
e*"' sin cos 







1 



Tpq{(f), a) is just a plane unitary transformation acting non-trivially only on states p and q. 

Let V G U{n). We want to find some T„^n-i(0, cr) such that T*^^_^V = V = [v^nxn, 
where = 0: 



n,n—l ' 



1 
1 
1 



o 



O 

COS0 

-e"^sin0 COS0 



^^11 

Vn-1,1 
Vnl 



Vl,n-1 Vin 

1^n—l,n—l 
^n.n— 1 ^nn 



SO 



n = Vn-l,n COS (j) + i;„„e siu (j). 



We consider all possibilites: 
Case 1: Vn-i,n = 0. Then we choose = 0, cr = 0, i.e., T'„_i,„(0, cr) = /„, and we obtain 



Case 2: Vn-i,n 0; '^nn = 0. Then choose = 7r/2, cr = 0. Obtain v'^_i ,^ — 0. 

Case 3: Vn-i,n + 0,^^^ 7^ 0. Write = r„-i,ne^^"-i'", 

-Qn-\,n + ^Inn and = tan~^ (-r„_i,„/r„„) . Obtain 



Choose a 



^n-i,n = COS • r„_i,„e 

^n— l,n 



+ sm • rnne 
+ tan ) r„„ cos ^e'^"-^'" = 0. 



Therefore, we have found T'n,n-i £ (n) such that 



n.n— 1 ' 



In 



<1 



n.n— 1 



n— 2,n 
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Similarly, we can find ^n,n-3, • • • , ^n,i such that 



ry-i* -IT- 

n,n—2 71,71—1 



an 



71— vi.n 



* 








71,71—1 nn 



-'- nl-'- nl ■ ■ ■ -'- n,n-2-'- n,n-l ^ 



Since W is unitary, we conclude f„i = Vn2 
a„ G M. Thus 



* 

Vnl 



* ■ 
: 

* 

^n,7i— 1 ^rtTi. 



w. 



-'nl-'n2 • • • -'- n,n-2-'- n.n-1^ 



Vn,n~i = and f„„ = e*"" = (i„ for some 










dn 



Now, applying the same technique to the remaining {n — 1) x {n — 1) undiagonalized ma- 
trix block (**) above, together with a simple induction argument, we obtain plane unitary 
transformation T„i, . . . , T„,„_i, T^-i^i, ... , Tn-i,n-2, ■■■ , T31, T32 and T21 such that 



rri^ rri^ rri^ rri^ rri^ rri^ rri^ rri^ t r 

-'21 -'si -'32 -'41 • • • -'71-1,1 • • • -'n-l,n-2-'7il • • • ^n,n-l^ 



di 



d2 




where dj = e*"^ for j = 1,2, 
Therefore 



n. 



^ — ^71,71-1; Tn^n-2 ■ ■ ■ ^nl^n-l,n-2 • • • ^n-1,1 • • • ^32^31^21-D 
n i—1 \ 



,1=1 j=i 



and (ITq) is proved. 



11 



References 



1. A. Barenco, C.H. Bennett, R. Cleve, D.P. Divincenzo, N. Margolus, P. Shor, T. Sleator, 
J.A. Smolin, and H. Weinfurter, Elementary gates for quantum computation, Physical 
Review A, 52(5) (1995), 3457-3467. 

2. R. Cleve, A. Ekert, C. Maccliiavello and M. Mosca, Quantum algorithms revisited, 
Proc. R. Soc. London A 454 (1998), 339-354. 

3. A. Ekcrt, Quantum interferometers as quantum computers, Phys. Scripta T76 (1998), 
218-222. 

4. A. Klappenecker, Computing with a quantum flavor, manuscript (2000). 

5. F.D. Murnaghan, The Unitary and Rotation Groups, Spartan, Washington, D.C., 1962. 

6. M. Reck, A. Zeilinger, H.J. Bernstein, and P. Bertani, Experimental realization of any 
discrete unitary operator. Physical Review Letters 73 (1994), 58-61. 



12 



